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First principles band-structure calculations for large supercells of Ba2Cu04-« and La2Cu04_a 
with different distributions and concentrations of oxygen vacancies show that the effective doping 
on copper sites strongly depends on where the vacancy is located. A vacancy within the Cu layer 
produces a weak doping effect while a vacancy located at an apical oxygen site acts as a stronger 
electron dopant on the copper layers and gradually brings the electronic structure close to that of 
La2-2,Sr 2! Cu04. These effects are robust and only depend marginally on lattice distortions. Our 
results show that deoxygenation can reduce the effect of traditional La/Sr or La/Nd substitutions. 
Our study clearly identifies location of the dopant in the crystal structure as an important factor in 
doping of the cuprate planes. 

PACS numbers: 74.72.-h, 71.15.Ap, 71. 15. Mb 



I. INTRODUCTION 

Recent work on high-T c cuprates has explored the role 
of oxygen defects and ordering of oxygen into super- 
structures. For instance, x-ray diffraction measurements 
on small patches of copper oxide superconductors reveal 
structures which may promote high temperature super- 
conductivity Other studies involving O vacancies in 
the high-T c cuprate Sr 2 Cu04_5 (SCO) have shown that 
the superconducting T c may become surprisingly large, 
reaching 95 K for certain heat treatments These su- 

perconductors have the same 214 structure as La2Cu04 
(LCO) and it is believed that the La substitution with 
Sr is associated with vacancy formation on the oxygen 
lattice. Interestingly, up to 30 percent of the oxygen 
might be missing in the superconducting samples [2j. 
This finding is at odds with the doping dependence of 
T c in La2- a; Sr a ;Cu04 (LSCO), since the maximum T c is 
found for a moderate hole doping x of about 0.15 holes 
per Cu. It is clear thus that it is important to under- 
stand the role of oxygen vacancies in doping the Cu-0 
planes, which are widely believed to be the seat of super- 
conductivity. So motivated, we examine in this article 
how location of the oxygen vacancies can affect the dop- 
ing of cuprate planes [4|, la] ■ An outline of this paper is as 
follows. In Sec. II, we present the details of the electronic 
structure computations in large supercells. The theoret- 
ical results are presented and discussed in Sec. Ill, and 
the conclusions arc summarized in Sec. IV. 



II. METHOD OF CALCULATIONS 

We specifically extract position dependent mechanisms 
by which the oxygen vacancies in Ba2Cu04 (BCO) and 
La2Cu04 modify the electronic structure. Since Ba is iso- 
electronic with Sr, our results for the valence electrons are 
also applicable to Sr 2 Cu04_5.[l, 0] Self-consistent first 
principles calculations are performed for (4,2,2) exten- 



sions of the basic cell [8] into supercells containing 16 
formula units with 112 atomic sites in total, as well as 
for smaller cells with 7 and 14 sites, by using the Linear 
Muffin-Tin Orbital (LMTO) method [|| and the Local 
Density Approximation (LDA) [l(| ■ While the LDA does 
not provide a satisfactory description of the electronic 
structure in underdoped cuprates because of inadequate 
treatment of correlations [Lj|, good agreement is found 
between LDA and various experiments [H-G]| for opti- 
mally and overdoped cuprates. The present calculations 
consider mostly the metallic regime where LDA is ex- 
pected to be a reasonable approximation. The converged 
self-consistent results are obtained using a mesh of 125 
k-points within the irreducible Brillouin zone. A precise 
tetrahedron method is used to determine the density-of- 
states (DOS) [l5j]. All sites in the cell are considered 
as non-equivalent throughout the self-consistency cycle. 
We have previously used supercell computations along 
these lines to search for weak ferromagnetism around 
Ba-clusters in La2- a; Ba a ;Cu04 (LBCO) p}. However, 
the present calculations are mostly non spin-polarized, 
since magnetism is not an issue here. In fact, earlier 
self-consistent spin-polarized band-structure calculations 
with 28-atom unit cells have revealed that oxygen vacan- 
cies strongly reduce the tendency toward antiferromag- 
netism in LCO [Hj]. The LMTO results for nonmagnetic 
La2Cu04 are in excellent agreement with full potential 
calculations [TtJ demonstrating the high quality of our 
basis set. 

Vacancies (V) on oxygen sites are introduced in such 
a way as to avoid clustering of vacancies. We consider 
configurations with number of vacancies ny equal to 1, 
5, 8, 9, 10 or 11. The corresponding oxygen deficiencies 
for Ba2Cu04_,5 are given by 5 = ny /16. Thus 5 ranges 
between 0.06 and 0.69. The maximum T c in the SCO sys- 
tem is reported to occur for 8 « 0.6, and another dome 
of high T c in an uninvestigated region of the phase dia- 
gram has been proposed [2| . The same type of supercells 
with ny equal to 5 and 9 are used in calculations with 
La replacing Ba in order to compare the trends for Cu-c? 
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band filling in BCO and LCO. The lattice constant is 
kept fixed in all calculations. The present computations 
are performed with the oxygen vacancies either in the 
planes or in apical sites. 



III. RESULTS AND DISCUSSION 

The Cu-d band filling can be monitored through in- 
spection of the DOS or the effective charge within the 
Cu muffin-tin spheres. The DOS at the Fermi energy, 
N(Ep), is dominated by Cu-d electrons in all high-T c 
cuprates, and the mechanism of superconductivity likely 
involves the Cu-3ci and 0-2p character on the Fermi sur- 
face (FS). Therefore, it is important to see to what ex- 
tent the d-band filling is affected by oxygen vacancies 
or other dopants. In the left insert to Fig. Q] is shown 
the total DOS for undoped BCO (for the same supercell 
as for doped cases), and with 9 vacancies either in the 
planes or at the apical positions. Despite the unavoid- 
able noise in the DOS, there is a clear d-band shift to 
lower energy (with respect the Fermi level) for vacancies 
on apical sites compared to the pure case (BCO), while 
the shift for vacancies in the planes is much smaller with 
even some small enhancement of the DOS above Ep. A 
DOS enhancement produced by O vacancies was found 
in coherent-potential-approximation (CPA) calculations 
[l8l Il9| where all the oxygen sites (planar and apical) had 
equivalent random occupation. Notably, rigid band like 
pictures [2(], [2l[ have often been invoked in describing the 
doping evolution of the overdoped and optimally doped 
cuprates. Our large supercell treatment here goes beyond 
the simple rigid band models or the CPA-type mean-field 
approaches [22j to elucidate the local electronic and mag- 
netic properties of the system. 

Insight into the underlying physics is provided by the 
main frame of Fig. [TJ which compares the partial DOS 
of planar and apical oxygen in the undoped parent com- 
pound. Note that since Ba donates one less electron than 
La, the Fermi level is shifted downward compared to the 
situation in LSCO. In this case, bands of both d x 2_ y 2 
and d z 2 character are present near the Fermi level, with 
the latter dominating the DOS as shown by right hand 
side insert in Fig.[TJ Consequently, the partial DOS of the 
apical O near the Fermi level is larger, since the electrons 
of these O atoms hybridize with the Cu d z 2 electrons. 

Figure [2] shows the total DOS for pristine BCO (for 
the same supercell as for doped cases), and with up to 9 
vacancies in the planes or up to 11 vacancies in the apical 
positions. Apical site deoxygenation shows a clear trend, 
shifting the d-band to lower energy. The effect is also con- 
sistent with the change of the number of valence electrons 
within the Cu muffin-tin spheres for different vacancy 
configurations, shown in Table HI The strikingly differ- 
ent DOS shapes of apical and planar O-p bands in BCO 
suggest that x-ray emission spectroscopy (XES) would 
be a good experimental method for detection of doping 
dependencies. Furthermore, the d band shift below Ep 
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FIG. 1: (Color online) Local oxygen-p part of the DOS on one 
apical and one planar oxygen site in pristine Ba32Cui60s4. 
Left insert: DOS near the Fermi energy for 3 configurations, 
Ba32Cui6064 (BCO) and vacancy doped Ba32Cui6055Vg, 
with the 9 oxygen vacancies in plane (P9) and apical (A9) po- 
sitions, respectively, all calculated from 125 k-points. Right 
insert: Local copper-d part of the DOS (total and d(z 2 ) con- 
tributions) near the Fermi energy. All the DOS plots include 
a small broadening of 5 meV. 



seen in Fig.[2]for the apical but not for planar O vacancies 
would be a key signature to look for with XPS. Notably, 
recent O K-edge x-ray absorption (XAS) experiments on 
SCO/LCO superlattices |23| have shown that the reduc- 
tion of the second hole peak in the unoccupied DOS of 
SCO provides a signature of apical oxygen vacancies in 
SCO. In fact, earlier XAS measurements on overdoped 
LSCO Hi] and electron energy-loss spectroscopy on SCO 
[25l ] have determined that such peaks belong to the api- 
cal oxygen unoccupied partial DOS. Clearly, changes in 
these unoccupied DOS peaks [HI mirror trends shown in 
Fig. [2] associated with the apical oxygen character in the 
occupied DOS. 

The Cu muffin-tin spheres have the same radius in 
all cases. For an increased number of apical vacancies 
in BCO there is a monotonic increase of the Cu charge 
so that about 0.1 electrons are added to the sphere for 
S r; 0.56. The fine details of the vacancy distribution 
do not seem to be important since results from two dif- 
ferent configurations (for 8 and 10 vacancies) give almost 
identical results. In BCO with vacancies within the Cu02 
planes, there is little modification of the Cu charge, about 
1/10 of the change due to apical vacancies. These small 
changes shown in Table Q] indicate less Cu charge for O 
vacancies within the Cu02 plane, i.e. there is an effective 
hole doping within the Cu band. 

The average valence charge of the oxygen vacancy Qv 
in Table I yields a simple physical picture of the differ- 
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TABLE I: Average valence charge of Cu (Qcu, within the 
atom sphere radius of 0.335ao) and its difference relative 
to the case without vacancy (AQcu) for different number 
of oxygen vacancies ny within planar and apical positions 



in Ba32Cui 6 6 



,V„ V (BCO) and La 32 Cui6064-n v V„ 



(LCO). The corresponding average valence charge of the oxy- 
gen vacancy (Qv within atomic sphere radii of 0.317ao and 
0.345ao for planar and apical sites, respectively) is also re- 
ported. Since there are in total 64 oxygen sites in the 112 
atom cell, the effective doping is S — nv/lG. Different distri- 
butions of 8 and 10 vacancies on apical sites for BCO show 
very little difference in the charge transfer. The cases nv = 1 
with or without lattice distortion are almost identical, see 
text. 



BCO BCO BCO LCO LCO LCO 

Qcu AQcu Qv Qcu AQcu Qv 



n v 



9 (plane) 
8 (plane) 
5 (plane) 
(plane) 



(apical) 
(apical) 
8 (apical) 

8 (apical) 

9 (apical) 

10 (apical) 

10 (apical) 

11 (apical) 



10.054 
10.055 
10.057 
10.062 
10.064 
10.071 
10.107 
10.142 
10.145 
10.164 
10.179 
10.180 
10.200 



-0.010 
-0.009 
-0.007 
-0.002 
0.000 
0.007 
0.043 
0.078 
0.081 
0.100 
0.115 
0.116 
0.137 



0.662 
0.655 
0.650 
0.640 

0.303 
0.313 
0.320 
0.308 
0.327 
0.326 
0.329 
0.336 



10.331 -0.063 0.889 

10.368 -0.026 0.877 

10.394 0.000 

10.450 0.056 0.623 

10.456 0.062 0.662 



ence between vacancies in planar and apical O-positions. 
When the oxygen vacancy is located in the Cu-0 planes, 
it acts as an electron trap (i.e. Qv is large), while when 
the vacancy is in the charge reservoir, the electrons are 
more attracted by the Cu-0 planes than by the vacancy. 
Hence, apical vacancies behave as expected from simple 
valence arguments: 2 ~ — > V + 2e, which dopes the 
planes with two electrons per vacancy. On the other 
hand, planar vacancies are anomalously charged, so : 
2 ~ — > V~ + e, resulting in a weaker electron doping. 
Note that in both cases, electrons are doped into the 
planes, but for planar O vacancies, they are more local- 
ized on the vacancy. Apical oxygen vacancies will thus 
mainly remove O-p bands in the interval where the apical 
O-DOS is high ranging from Ep down to 1.5 eV below 
Ep. In contrast, the DOS reduction due to planar oxy- 
gen vacancies will be spread over a wider energy range 
because of the wider partial DOS for planar oxygen sites 
shown by Fig. [TJ 

We also carried out computations on LCO to further 
demonstrate the robustness of our results. In particular 
as illustrated in Table HI the electron doping from apical 
O- vacancies in LCO is comparable to that of BCO. More- 
over, to estimate the effect of lattice relaxation around 
vacancies, calculations were made with a vacancy on one 
apical oxygen position for BCO, where a large lattice 
relaxation is imposed. The two atoms surrounding the 
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FIG. 2: (Color online) DOS near the Fermi energy for var- 
ious vacancy configurations in Ba32Cui6064-n- BCO is the 
undoped case, "Pn" refers to n vacancies in the plane and 
"An" refers to n apical vacancies. The DOS plots include a 
small broadening of 5 meV. The Fermi energy is taken as zero 
in all cases. Black dashed lines show the DOS for different ar- 
rangements of the same number and type of O vacancies. This 
illustrates the robustness of the main trends. When n > 8 for 
the apical vacancies, the arrows indicate the onset of the sharp 
rise associated with the d,i band below the Fermi level. 



missing O site, namely the Ba above (along the c axis) 
and the Cu below were allowed to fill in the empty space 
by moving 0.03ao towards the center of the missing O. 
The amplitude of these distortions is rather large, corre- 
sponding to a thermal distortion amplitude for a temper- 
ature of several hundred degrees Kelvin [26| . This distor- 
tion is found to make very little change in the number of 
electrons on the Cu sites. No Cu atom changes its charge 
more than 0.005 electrons. However, the Cu atom next 
to the vacancy gains 0.05 electrons compared to Cu in 
ideal BCO, and this is independent of relaxation. The 
Ba atom looses 0.08 electrons and the vacancy (with the 
same volume as the O atom) gains 0.10 electrons. Other 
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sites typically change their charge by ±0.01 electrons or 
less. Therefore, our finding of net charge transfer due to 
high vacancy concentration is essentially independent of 
lattice relaxations. 

To better understand the role of apical vacancies in 
doping the planes, we also present some calculations for 
a related compound with even more apical O vacancies: 
Ba2CuC>3 (BCO-3) where one layer of apicals is miss- 
ing [13. Figure |3] compares the DOS of BCO-3 to sto- 
ichiometric Ba 2 Cu0 4 (BCO-4) and La 2 Cu0 4 (LCO-4), 
demonstrating that the DOS of BCO-3 and LCO-4 are 
almost identical near Ep - i.e., for the d x 2_ y 2 band. At 
lower energies, from 1 eV below Ep, there arc differences 
because of the removal of three hybridized O-p bands. In 
contrast, the DOS of BCO-4 is very different even near 
Ep. These DOS functions are consistent with the trend 
shown in Fig. [2] despite the different compositions and 
orderings. Disorder has some effect on the details but not 
on the main features of the DOS. The Cu valence charge 
in BCO-3, 10.25 electrons, is consistent with the results 
in Table U for an extrapolated concentration of 16 apical 
vacancies. Thus, the DOS and FS of the relevant Cu- 
O single band of BCO-3 and LCO-4 behave in the same 
way, but for an effective Cu charge that is slightly smaller 
in BCO-3: about 10.25 vs. 10.39 in LCO-4, see Table Q] 
There is a striking similarity of the FS of BCO-3 (shown 
in the upper inset of Fig. [3]) and the well-known FS of 
LCO-4 (see ref. [28j]). Furthermore, it can be verified 
that the FS's of the two systems have an almost iden- 
tical evolution with the number of holes p in the single 
band at Ep as shown in the insert of Fig. [3] for BCO-3. 
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FIG. 3: (Color online) DOS near the Fermi energy for 
Ba2Cu03, Ba2Cu04 and La2Cu04. Insert (a) shows the evo- 
lution of the Fermi surface in the k z = 0- plane in Ba2Cu03 as 
a function of the rigid-band doping for 0, 0.15 and 0.30 holes 
per unit cell, while insert (b) shows how the Fermi level shifts 
above the d z 2 band onset as vacancies are added. 

We can now understand that the role of apical vacancy 



doping is to push the Fermi level from the middle of the 
d z 2 band for xiy = to the middle of the d x 2_ y 2 band for 
ny = 16. For riy > 8, the sharp rise associated with the 
d z 2 band is located below the Fermi level as shown by 
the arrows in Fig. [51 We plot the energy shift between 
the d z 2 onset and Ep as a function of S = ny/16 in 
insert (b) of Fig. |3] Remarkably, the shift for BCO-3 
also fits on the same line. A linear shift is expected if the 
DOS is approximately constant, as in Fig. [3] Hence, we 
find that in BCO for S between 0.5 and 1, the vacancies 
dope electrons into an essentially empty d x 2_ y 2 band, and 
the physics should be similar to that of strongly hole- 
overdoped LCO. 

In comparing our results to experiments on the 
Sr2Cu04_5 system, we note that Gao et al. [|| conclude 
that the missing oxygens are outside the CuO plane and 
that the structure remains stable even at this high dop- 
ing [2^|. These conclusion are also corroborated by our 
calculated total energies for having 1 apical or 1 planar 
vacancy in the undistorted 112 site supercell, since the 
total energy is considerably lower for the apical position, 
of the order of 1.5 eV. Vacancy ordering leading to long 
range periodicity is suspected to be important for boost- 
ing T c to high values [l], [30l - [32] | . Orderings of apicals or 
other defects, leading to long-range periodicity, are ex- 
pected to induce a lowering of the total energy because 
of the appearance of a partial gap at Ep [26j. 

Since each apical vacancy contributes two electrons to 
the Cu02 plane, the addition of apical vacancies could 
lead to a new route for electron-doping LCO. We note 
that O vacancies in the reservoir layers also provide elec- 
tron doping in YBa2Cu307_ y (YBCO) |33| . even though 
YBCO is structurally very different from both LCO and 
BCO. The doping for high-T c s in BCO cuprates corre- 
sponds to a nearly empty d x 2_ y 2 band, consistent with 
the suggestion of a second superconducting dome. It is 
interesting to note that dopings for high-T c in pnictides 
often seem to be associated with the Van Hove singular- 
ity at a band onset [34j], although the reason for this is 
not clear. In the case of BCO, the possibility of altered 
d z 2-d a .2_ J( 2-orbital character near Ep due to vacancy or- 
dering could also modify the superconducting properties. 



IV. CONCLUSION 

Our study shows that properties of copper oxides 
high temperature superconductors depend dramatically 
on the concentration and position of oxygen vacancies. 
Our results indicate that manipulation of apical oxygen 
vacancies is an effective pathway for optimizing the elec- 
tronic structure for higher T c . The similarity between 
the band structure for apical vacancy doped BCO and 
the commonly doped LCO system is striking, but con- 
sistent with the notion that the high T c in Sr2Cu04_^ is 
due to another " dome" of superconductivity with a dif- 
ferent hole doping on Cu-sites. The slight difference in 
Cu valence charge between the two systems can affect 
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the magnetic response despite a very similar Fermi sur- 
face. Another implication of our study is that O vacancy 
formation can compensate the doping resulting through 
conventional La substitutions. 
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